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ABSTRACT 

Older and more recent pieces of observational evidence suggest a strong connection 
between QSOs and galaxies; in particular, the recently discovered correlation between 
black hole and galactic bulge masses suggests that QSO activity is directly connected 
to the formation of galactic bulges. The cosmological problem of QSO formation is 
analyzed in the framework of an analytical model for galaxy formation; for the first 
time a joint comparison with galaxy and QSO observables is performed. In this model 
it is assumed that the same physical variable which determines galaxy morphology is 
able to modulate the mass of the black hole responsible for QSO activity. Both halo 
spin and the occurence of a major merger are considered as candidates to this role. The 
predictions of the model are compared to available data for the type-dependent galaxy 
mass functions, the star-formation history of elliptical galaxies, the QSO luminosity 
function and its evolution (including the obscured objects contributing to the hard- 
X-ray background), the mass function of dormant black holes and the distribution of 
black-hole - bulge mass ratios. A good agreement with observations is obtained if the 
halo spin modulates the efficiency of black-hole formation, and if the galactic halos at 
z = have shone in an inverted order with respect to the hierarchical one (i.e., stars 
and black holes in bigger galactic halos have formed before those in smaller ones) . This 
inversion of hierarchical order for galaxy formation, which reconciles galaxy formation 
with QSO evolution, is consistent with many pieces of observational evidence. 

Key words: cosmology: theory - galaxies: formation - quasars: general - dark matter 
- large-scale structure of Universe 



1 INTRODUCTION 

High-redshift quasi-stellar objects (QSOs) have been for 
a long time the only probe to the high-redshift Universe 
at 1 ^ 2 4. However, their potential power in con- 
straining cosmological models has always been hampered 
by their complexity as a physical phenomenon: they are 
thought to be powered by huge black holes (BHs), of mass 



10" 



10 Mq, hosted at the center of proto-galactic ha- 



los, so that their activity couples very different scales, from 
fractions of pc to Mpc. The problem is then unsuitable for 
numerical studies, and the modeling of QSOs in a cosmolog- 
ical framework must rely on analytical approximations, typ- 
ically based on uncertain hypotheses or poorly constrained 
parameters. 

The QSO population is thought to be made of many 
generations of short-lived events, with lifetimes ranging from 
a few lO'^ to 10* yr, as longer lifetimes would imply few 
very large dormant BHs, with masses ~ 10^^ M©, which is 
contrary to the observational evidence (see, e.g., Cavaliere & 



Padovani 1988) . This fact has two implications, which reveal 
the deep connection between QSOs and galaxies: firstly, the 
number density of expected dormant objects matches that 
of bright galaxies; secondly, to evolve on cosmological time 
scales as observed, the various generations of QSOs must be 
coordinated by a process {great coordinator), which is likely 
to be that of galaxy formation. 

In fact, while QSOs have always been addressed as a 
separate field in cosmology, both theoretically and observa- 
tionally, their connection with galaxies and galaxy formation 
is receiving ever growing evidence from low and high redshift 
observations. There is no evidence of QSO activity outside 
galaxies, and the direct observation of high-redshift galax- 
ies has revealed that large spheroidal galaxies are the most 
common hosts for bright QSOs (see, e.g., Hall & Green 1998; 
McLure et al. 1998). Massive dark objects, interpretable as 
large dormant BHs, are routinely found in nearby spheroids 
(Kormendy & Richstone 1995; Ford et al. 1997; Magorrian 
et al. 1998; van der Marel 1998; Ho 1998; Wandel 1998). The 
estimates of their masses are still affected by many systemat- 
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ics, but two resulting evidences seem robust: (i) the mass of 
the massive dark object is correlated to the mass of the bulge 
component, (ii) this correlation has a scatter of about one 
order of magnitude. Finally, the light history of QSOs has 
an interesting resemblance with the star-formation history 
of galaxies (Boyle & Terlevich 1998; Cavaliere & Vittorini 
1998). 

The presence of these large dormant BHs in the cores 
of the bulges of nearby galaxies is a key quantity for test- 
ing the BH paradigm of QSOs (Soltan 1982). A previous 
paper (Salucci et al. 1998a; hereafter paper I) has been ded- 
icated to finding the mass function of dormant BHs, us- 
ing up-to-date observations and including the contribution 
of a population of heavily obscured objects (type II AGN) 
which are revealed by their contribution to the hard- X-ray 
background (see also Iwasawa & Fabian 1998). The mass 
function of dormant BHs in nearby galaxies has been es- 
timated using two methods. The first exploits the recently 
discovered correlation of BH and bulge masses, and consists 
in convolving the mass function of the galactic bulges with 
a BH-buIge relation inferred from observation. The second 
method relies on the correlation between radio power and 
BH mass (see also Franceschini, Vercellone & Fabian 1998). 
These two estimates agree with the mass function of the 
accreted matter inferred from the QSO emission (including 
obscured objects), giving direct support to the QSO galaxy 
connection. The local mass density in BHs turnes out to be 
~ 6.5 10^ Mq Mpc"^ for Ho = 70 km/s/Mpc 

The results found in paper I reveal a dichotomy be- 
tween large BHs {M > 10* Mq) and small BHs. The former 
are hosted in elliptical galaxies, shine only once near to the 
Eddington rate, and tend to be not obscured, while the lat- 
ter are found in the bulges of spiral galaxies, shine with 
a lower efficiency, can be reactivated by interactions and 
are frequently obscured. In Salucci et al. (1998b) the mass 
function of BHs in spirals has been constrained through the 
use of several hundred rotation curves of spirals; no BH is 
detected, and the upper limits thus obtained constrain the 
numerous late spirals to host a negligible amount of mass in 
BHs. 

Within the framework of hierarchical cosmological mod- 
els, it is possible to make prediction on the number and 
properties of dark-matter (DM) halos. Besides, to predict 
the statistical properties of QSOs it is necessary to make as- 
sumptions on the probability that a BH forms inside a DM 
halo, and the efficiency with which such a BH radiates en- 
ergy. Arguments in favor of BH formation in normal galactic 
halos are given, for instance, in Rees (1984) and Haehnelt & 
Rees (1993). Efstathiou & Rees (1988) made the assumption 
of constant ratio between halo and BH mass and of radia- 
tion at the Eddington limit, to conclude that the standard 
cold dark matter (CDM) model could reproduce the QSO 
luminosity function. A similar but more refined procedure 
was used by Haehnelt & Rees (1993): they assumed that 
BHs form with an efficiency which increases with central 
halo density and halo virial velocity, so as the QSO activ- 
ity reaches a maximum at z ~ 3 as observed. Haehnelt, 
Natarajan & Rees (1997) considered the case in which most 
mass is accreted during the quiescent phase, while Catta- 
neo, Haehnelt & Rees (1999) analyzed the implications of 
the BH-bulge correlation, especially on galaxy mergings. A 
related approach was used by Katz et al. (1994), who used 



N-body simulations to check whether the CDM model could 
give a sufficient number of suitable halos to justify QSOs 
at z <: 3. Carlberg (1990), and more recently Krivitsky & 
Kontorovich (1998), estimated the number of QSOs by as- 
suming them to be related to galaxy mergings. Predictions 
on the cosmological evolution of QSOs were given by Haiman 
& Menou (1998) and Percival & Miller (1999). Eisenstein & 
Loeb (1995a,b) followed in some detail the dissipation of an- 
gular momentum of gas infalling inside a DM halo, modeling 
the collapsing halos as homogeneous ellipsoids and seeking 
for suitable conditions for a BH to form: they concluded that 
seed BHs (with mass ^ lO'' Mq) can form at rather large 
redshifts (z > 5), giving rise to QSO activity when included 
in a protogalaxy (see also Loeb 1993). Cavaliere & Vittorini 
(1998) argued that QSOs at z <: 3 are connected to the 
formation of new galactic halos, while the already- formed 
BHs are reactivated by galaxy interactions in newly-formed 
groups, a mechanism which is effective at z ^ 3 (see, e.g., 
Monaco et al. 1994, and references therein). 

This paper aims to construct, for the first time, an ana- 
lytical model for joint QSO and galaxy formation, which re- 
produces the main observables relative to both populations. 
Mergers and/or halo spin are proposed as possible physical 
variables responsible for both galactic morphology and BH 
formation. The analytical model is designed to reproduce 
the mass function of galactic halos for different broad mor- 
phological classes, the star formation history of ellipticals, 
the QSO luminosity function and its redshift evolution, the 
mass function of dormant BHs and the BH-bulge relation. 
The plan of the paper is as follows: Section 2 presents the 
analytical model for galaxy formation, based on the joint 
distribution of halo mass, spin or last merger, and forma- 
tion time. In Section 3 the model is shown to reproduce the 
estimated halo mass function of galaxies, divided into broad 
morphological classes, and the star-formation history of el- 
lipticals. In Section 4 the model is compared to the observed 
QSO luminosity function. Section 5 addresses the dormant 
BH masses and their relation to bulge masses. Section 6 
contains a summary and some final remarks. 

A Hubble constant of Ho = 50 km Mpc~^ will be 
used when discussing Einstein-de Sitter (EdS) models, while 
a value of 70 will be used for the open model or the model 
with cosmological constant. 



2 AN ANALYTICAL MODEL FOR GALAXY 
AND QSO FORMATION 

Galaxy formation in hierarchical models is usually addressed 
by means of semi-analytical techniques, in which the abun- 
dance of DM halos and their merging histories are inferred 
either from N-body simulations or from the extended PS for- 
malism (Bond et al. 1991; Bower 1991; Lacey & Cole 1993), 
and the history of gas in halos is described through a set of 
simplified rules for gas cooling, star formation, feedback pro- 
cesses, galaxy mergings etc. (see, e.g. White 1993; Cole et al. 
1994). Noteworthy, the main qualitative conclusions of these 
semi-analytic models can be reached by means of simple an- 
alytic arguments based on the PS mass function (White & 
Rees 1978; White & Frenk 1991). It is then reasonable, at 
this stage, to construct a simple analytical model for joint 
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galaxy and QSO formation, leaving detailed calculations to 
further analysis. 

The power spectra considered in the present paper are 
the CDM-like ones described by the parameterization of Efs- 
tathiou, Bond & White (1993), with the shape parameter de- 
fined as r = Qh, where Q is the cosmological density param- 
eter and h = //o/(100 km/s/Mpc). Three models have been 
analyzed, namely an EdS model with — 1 and h — 0.5 
(so that r — 0.5), a flat low-density model with fl = 0.3 
and cosmological constant {h = 0.7 and F = 0.21), and an 
open model with D. = 0.3 and h = 0.7 (F = 0.21). Follow- 
ing Eke et al. (1998), the normalization has been fixed so 
that the standard deviation of the initial density (linearly 
extrapolated to z = 0) on a 8 Mpc sphere, as, is 0.7 
for the EdS model, and 1 in the other cases. For the sake 
of brevity, the results of the open model, which are similar 
to the cosmological constant case (hereafter called Lambda 
model), are not shown. 



2.1 Dark-matter halos in hierarchical Universes 

The mass function of DM halos at a given redshift is repro- 
duced by the well-known PS formula (see Monaco 1998 for 
a recent review on the mass function) : 



nps{MH;z)dMH = 

Po 1 
Mh ^/2^ 



(1) 



exp 



5l/b{zf 
2A 



dA 



dMn, 



where Mh is the halo mass, po is the background density, 
A(Mif) = ((5^) is the mass variance at the scale correspond- 
ing to Mh (it is usually denoted as ct^), and 5^ is a threshold 
parameter. Following Monaco (1998, 1999) and Governato 
et al (1998), the Sc parameter is set to 1.5 if the cosmolog- 
ical density is f2 = 1, and 1.69 if f2 < 1 (with or without 
cosmological constant). The time-dependent function b{z) is 
the linear growing mode, normalized as b{z = 0) = 1. This 
quantity is related to the scale factor a{z) = 1/(1 + z) as 
follows: b{z) oc a{z) at high redshift, and b{z) — a{z) in the 
n = 1 cosmology. The expressions for b{z) for cosmologies 
different from EdS are given, for instance, in Monaco (1998). 
The critical mass M, (z) at a given redshift is defined as the 
mass at which the variance (linearly extrapolated to z) is 
equal to A{M4z)) = {Sc/b{z)f. 

Galaxies do not form in every DM halo: large halos have 
large cooling times, which can become larger than the dy- 
namical time, in which case baryons cannot gather into a sin- 
gle unit. Following White & Rees (1978) and White (1993), 
the virial temperature Th of a halo of mass Mh and density 
Qh is assumed to scale as Th oc M'^^ g^J'^ , while the cooling 
time depends on the efficiency of cooling Acooiing(r_ff) in the 
following way: tcooi oc Th Qh^ Acooiing{TH)~^ ■ In the relevant 
range of temperature, from ~ 10* to ~ 10^'^ K, the cooling 
efficiency scales roughly as Acooiing(2if ) oc T^^^^; then the 

— 1/2 

cooling time scales as fcooi oc MhQh ■ Finally, the dynami- 

— 1/2 

cal time of the halo scales as tdyn cc Qh ■ Then, a condition 
^cooi < idyn translates into a redshift-independent cutoff of 
the mass function. 

The largest halo mass in which gas can effectively cool 
down is called M^ooi, and left as a free parameter. Assuming 
that galaxies obey a relation Mh oc , with P ~ 0.5 — 1 



(see Section 3.1; /3 depends on the morphological type, we 
use the value relative to elliptical galaxies), we model the 
astrophysical cutoff so as to reproduce the observed cutoff 
in the galaxy luminosity function: 



Ccooi(A/h) = exp(-(M^f/Afeooi)'/''). 



(2) 



This astrophysical cutoff is not the only feature that 
distinguishes the assembly of galaxies from that of DM: 
the merging of baryonic matter is regulated both by cool- 
ing and feedback, when gaseous clumps merge, and by non- 
dissipative merging of clumps made of stars. To address the 
formation of a galaxy, it is necessary to distinguish between 
early-forming halos, whose density is high enough for sub- 
clumps to merge into a single galaxy, and late- forming halos, 
which host small groups of galaxies (see Section 2.4). We 
consider also a further parameter p, responsible for mor- 
phology. The number density uh of halos with mass Mh, 
morphological parameter p and formation redshift can be 
expressed with great generality as follows: 

nH{MH,p,zf)dMHdpdzf = nps{MH\zf) x Ccooi(MH)dAfe x 
Pp{p\MH,Zf)dp X Pf{zf\MH)dzf. (3) 

The following subsections are dedicated to finding suit- 
able approximations for the distributions Pp{p\AdH, Zf), 
with p equal to the halo spin A or merging fraction /, and 

Pf{zf\MH). 



2.2 The joint mass-spin function 

Spin is acquired by a halo during the mildly non linear 
regime, when the tidal coupling between the (proto-)halo 
and the external mass distribution is effective. After de- 
coupling, the halo inertia moment becomes negligible, and 
spin stops growing, at least as long as the structure remains 
isolated. It is then possible to give analytical estimates of 
halo spins (Peebles 1969, White 1984, Heavens & Peacock 
1988, Eisenstein & Loeb 1995a, Catelan & Theuns 1996). 
The spin distribution of halos is also determined by means 
of N-body simulations (Barnes & Efstathiou 1988; Zurek, 
Quinn & Salmon 1988; Warren et al. 1992; Ueda et al. 1994; 
Lemson & Kauffman 1998). The comparison of analytical 
and numerical estimates leads to the following conclusions: 
(i) the spin distributions given by different authors and with 
different methods are approximately consistent with each 
other; (ii) the functional dependences of spin on other halo 
parameters, predicted by means of analytical arguments, are 
consistent with N-body results; (iii) the spin distribution is 
wide, and nearly lognormal in shape. 

The dimensionless spin parameter A is defined as: 



A = LE^''^G~^M~^''^, 



(4) 



where L is the final angular momentum of the halo, E is 
its total energy, G is the gravitational constant and Mh is 
the halo mass. The spin parameter A is nearly independent 
of everything, except a weak dependence of its mean value 
/ = (log A) on halo mass: 



1(Mh/M,{z)) = lo~OLx \og{MH/M,{z)). 



(5) 



With = log(0.04), the average value of the spin parame- 
ter, for a set of halos with mass not much smaller than Af, , 
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is about 0.05, the value given by N-body simulations. The 
exponent a\ is in the range 0.1 - 0.2; the value 0.15 will be 
used in the following. This trend is theoretically explained 
by a dependence of angular momentum on the initial height 
of the peak from which the structure comes from (Catelan 
& Theuns 1996), and has been revealed in the N-body sim- 
ulations of Ueda et al. (1994) and Cole & Lacey (1996). 
The minus sign implies that rare, massive halos (relative to 
Mt(z)) tend to have a lower spin. 

The PDF of the A parameter, as given by Barnes & 
Efstathiou (1988), Warren et al. (1992), Ueda et al. (1994), 
Eisenstein & Loeb (1995), Catelan & Theuns (1996) and 
Lemson & Kauffman (1997), is well approximated by the 
following lognormal distribution: 



P4X\MH,z)d\ = 



exp 



2^1 



dlogA, (6) 



where a\ — 0.3 and I is given by Eq. |^. Note that some 
authors, as Mo, Mao & White (1998), use for a\ a value of 
0.21, which is 30% smaller than the best-fit one used here. 



2.3 The joint mass-merging function 

Merging histories of DM halos are correctly reproduced by 
means of an extension of the PS approach to the mass func- 
tion problem (Bond et al. 1991; Bower 1991; Lacey & Cole 
1993). The probability that a halo of mass Mi (correspond- 
ing to a variance Ai) is included into a halo of mass M2 
(corresponding to a variance A2) after a time interval A log 6, 
where the growing mode hiz) is used as time variable, is: 



P(Afi ^ Ma; z)dMiA log b ■- 



2n b(z 



.(Ai-A2)-''''x 



exp 



(Alogb)^ 
b{z) I 2(Ai-A2) 



dAiAlogfa. 



(7) 



Denoting the merging fraction / as the ratio between the 
masses of the merging and the final clumps, / = M1/M2, 
it is straightforward from Eq. |^ to obtain the probability 
Pf{f;MH,z) that a halo (with final mass Mr = M2) at 
redshift z has experienced a merging, within a time interval 
A log 6, with another halo of mass (M2 — Mi). In the fol- 
lowing. Ml will denote the major progenitor of the halo, so 
that Ml > M2/2 and 1/2 > / > 1. The timescale Alogfo is 
left as a free parameter. 

The function Pf{f; Mr, z) is straightforward to obtain, 
but its expression is cumbersome; it is shown in Fig. 1 in the 
case Mh = M* and A log 6 = 0.1. For / -> 0.5 the proba- 
bility does not vanish but saturates to a finite value, with a 
very fiat shape, which is not unexpected as the merging of 
two comparable masses is not a rare event. In Section 2.6 it 
will be shown that this is a problem for using the merging 
fraction as a variable which modulates the efficiency of BH 
formation. 

In the simple case of power-law power spectra, P{k) oc 
fc", the probability of having had a major merging, / < 
/c (where /c is a threshold) depends on Mh / Mt,{z) in the 
following way: 



2 - 



1 - 



0.4 



o.e 



8 



Figure 1. The Pf{f) distribution derived from Eq. |^, for Mjj = 
M* and A log b = 0.1. The merging fraction / is the ratio between 
the masses of the two merging clumps (larger over smaller). 



/ \("+3)/6 

P(/ < ,fc;A/H,^) OC I I X 



\M,{z 



(8) 



exp 



-const X 



M-M 



(n+3)/3 



X A log b 



Eq. ^ shows that the probability of having experienced a 
mayor merging event grows with mass, in a way similar to 
that of the spin parameter: the exponent is small again, of or- 
der 0.3 if the spectral index n is about —2. At variance with 
the spin case, this trend is not univocal: very large halos, 
with masses Mh S> M, do not experience major mergings 
(they are just rare!); this fact does not have a great im- 
portance, as these very rare objects turn out to be mostly 
irrelevant. 

2.4 The redshift of halo formation and shining 

The assembly of baryonic structures like galaxies does not 
follow the assembly of DM halos: a DM halo can host more 
than one galaxy, and galaxies inside DM halos can subse- 
quently merge. At the present time, most galaxies are not 
isolated structures, but are contained in larger halos, such 
as groups and clusters. As a consequence, most present-day 
galactic halos are not related to the isolated halos described 
by the PS mass function; those halos have formed, and were 
present "in the PS sense" at high redshift, and have sub- 
sequently gathered into larger structures. The cores of such 
halos have retained their identity as they had much higher 
densities than the groups in which they have fallen. 

In order to describe the formation of such galactic ha- 
los, it is necessary to have some information on their dy- 
namical history. This information is not contained in the 
simple PS mass function, but can be obtained in the ex- 
tended PS formalism, for instance with the semi-analytical 
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merging trees technique (Lacey & Cole 1993), in which the 
merging histories of a sample of halos are given. A simpler, 
analytical procedure was recently proposed by Percival & 
Miller (1999), who tested its validity against large N-body 
simulations. The extended PS theory gives the probability 
P{Mh\z) that a halo of mass Mh is present at a given red- 
shift z; it is possible to invert such probability to obtain the 
probability Pf{zf\MH) that a halo of mass Mh forms at the 
redshift zf. 

PAzAMH)dz, = -^exp j --dz,. (9) 

As in Section 2.3, the variance A{Mh /M,(z)) yields a 
dependence on the halo mass in units of the critical mass 
M*(2). As the curve is peaked on A ~ 1, objects form pref- 
erentially when their mass is not very different from M« , and 
then larger masses form at larger times. Another characteris- 
tic of Eq. ^ is that a large number of halos of mass ~ 
are predicted to form at small redshift, Zf < 1. Such ob- 
jects can be either galaxies or small groups, depending on 
whether the baryonic substructures manage to merge into a 
single entity. The cross section for dissipationless merging is 
proportional to the square of the halo density, which scales 
as the cosmological density. As a consequence, the proba- 
bility that a halo is going to host a single galaxy (at z=0) 
is suppressed at lower redshift as {1 + zj)^. Then, in order 
to pick up the galactic halos, the Percival & Miller (1999) 
probability of formation redshift is multiplied by: 

('+(^)y'^ ™ 

The reference redshift zo is left as a free parameter. It is 
useful to stress that this condition, related to dissipationless 
merging, is a complement to the cooling condition given in 
Eq|. 

The formation of stars in spheroids and the bright phase 
of QSOs are assumed to be close in time. The high metallic- 
ities inferred from QSO spectra suggest that significant star 
formation in the host galaxies has already taken place before 
the bright QSO phase (Hamann & Ferland 1993). Besides, 
QSO activity cannot easily have place much later than the 
formation of bulge stars, when cold gas is almost absent and 
the feeding of BHs is difficult. 

QSO activity does not take place immediately at the dy- 
namical formation time of a generic (galactic) DM halo. We 
assume that the 'shining phase of the halo, i.e. when QSO 
activity takes place and star formation is already turned on, 
is delayed by a time tdciay with respect to the dynamical 
formation of the halo. This delay is assumed to be longer 
for smaller halos, such as to invert the hierarchical order for 
halo shining. In this case, the brighter QSOs can shine be- 
fore the fainter ones, as observed. An alternative scenario 
able to reproduce the QSO luminosity evolution requires 
lower efficiency of BH formation with increasing time and, 
as a consequence, in larger halos (Haehnelt & Rees 1993). 
This would induce an anti-correlation between bulge and BH 
mass, contrary to the observational evidence. 

* In the present paper the phrases 'shining of the halo' and 
'galaxy formation' are used as synonyms. 



As mentioned above, metallicity studies of the QSO en- 
vironments show that significant star formation in the host 
galaxies has occurred before the QSO shining phase. Thus 
we can infer that star formation in larger early-type galaxies 
is turned on more rapidly. In other words, the QSO evolu- 
tion marks the history of the star formation rate in early 
type galaxies (Silva et al. 1999). Shorter time scales for star 
formation in massive ellipticals have been suggested by Mat- 
teucci (1994) and Bressan, Chiosi & Tantalo (1996), on the 
basis of the chemical evolution of star populations. 

The shorter delay of the QSO shining phase for larger 
BHs may be ascribed to several mechanisms. For instance, 
more powerful objects can be able to remove earlier, and 
from a larger solid angle, the dust surrounding the circum- 
nuclear regions. The gas in larger galaxies may gather more 
easily in the core of the host halo, because of lower angu- 
lar momentum, and cool down, thus accreting on a seed 
BH. Smaller halos may accrete gas later in secondary infall 
events. An additional possibility is that large halos are en- 
dowed of dense peaks which may merge rapidly to produce 
supermassive BHs. 

Since the details of processes involving baryons are not 
fully understood, we presently explore the mass-dependent 
delay hypothesis as a heuristic (and parametric) guess, 
which helps to reconcile the 'anti-hierarchical' evolution of 
QSOs with the bulge-BH relation, and is consistent with 
evidences coming from the study of stellar populations in 
ellipticals. 

If t/ = t{zf) is the dynamical formation time, the "shin- 
ing" time of the halo is tah = tf + tdeiay , and its "shining" 
redshift z^h ~ z{tsi^). The halo mass Mh will denote in the 
following the mass of the halo at the shining time tsh', the 
number density of halos will be calculated at the shining 
redshift. We have found it convenient to parameterize the 
dependence of the delay time on the halo mass as follows: 

tdolay(MH) = log (iQ^ f ~ f -log M§ J _^ ^^tf^ ^ ^^^^ 

Here t/ is the delay of a galaxy of halo mass larger than 
M§^, corresponding to an L« elliptical (see Section 3.1). 
The parameterization is such that while large halos are 
delayed hy tf, smaller ones are delayed proportionally to 
log(MH/Mf J. Eq. ^ multiplied by the cutoff given by 
Eq. |l^, is evaluated at the shining redshift Zsh ~ z{tBh)', 
in other words, the curve is shifted in time by an amount 

^dclay • 

Notably, the final number of objects is calculated by in- 
tegrating the contributions at various redshifts. This is not 
strictly correct, as small halos at a time can be part of larger 
halos at a following time; this introduces an uncertainty in 
the normalization of the number of halos. However, the effect 
is likely to be modest, especially if large halos shine before 
smaller ones. The use of semi-analytical techniques, based 
on the merging trees, would solve this problem, and would 
allow a more detailed description of the delay time. More- 
over, it would allow to relax the assumption, implicit in this 
approach, that the merging of the halos of already-formed 
spheroids is negligible. As mentioned at the beginning of 
this Section, this analytical approach is supposed to catch 
the most important dependences, while further refinements 
are left to future work. 

To quantify and visualize the inversion of hierarchical 
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Log (MJ Log M„ (Mg) 



Figure 2. Left panels: redshift of halo shining for elliptical galaxies, as a function of the halo mass. The lookback time in Gyr is 
given on the right axis. The thick continuous line gives the redshift at which 50% of halos have shone; the thick dashed lines are what is 
obtained if there is no delay. The dot-dashed, dashed and dotted lines give the redshift interval within which 68%, 95% and 99% of halos 
shine. The shaded areas highlight those galactic halos whose abundance is very small (their mass is larger than 3 times the halo mass 
corresponding to an L» elliptical). Right panels: delay time tdclay a function of halo mass Mh- The Eds (upper panels) and Lambda 
(lower panels) models are shown. 



order, it is useful to consider how the abundance of galactic 
halos of fixed mass grows in time. This is not the same as 
plotting Eq. ^ which does not take into account the different 
abundance of halos of given mass at different times. The 
abundance of galactic halos is presented below (Eq. ^2|) . The 
left panels in Fig. 2 show the distribution of shining redshifts 
for halos of fixed mass which host ellipticals (with spin 
threshold) at z — 0, for the EdS and Lambda cosmologies. 
The thick line shows the redshift at which half of the z = 
halos have shined, the other lines show the redshift intervals 



within which 68%, 95% and 99% of halos shine. The right 
panels show the delay time fdeiay(A/H) used. The best-fit 
parameters have been used, see Table 2 and Section 2.7. 
For comparison, the thick dashed line shows the 50% line 
obtained by assuming no delay. The inversion of hierarchical 
order is visible in the change of slope of the "iso-shining" 
lines at moderate and small halo masses. 

Fig. 2. predicts that many big ellipticals form at high 
redshift; because of bias, these will preferentially end up in 
clusters. This is consistent with the observational evidence. 
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based on the colour-magnitude relation (Bower, Lucey & 
Ellis 1992; Ellis et al. 1997; Kodama et al. 1998; but see 
also Shioya & Bekki 1998), the tightness of the fundamental 
plane (Renzini & Ciotti 1993; van Dokkum et al. 1998), and 
the Mgi, — (Jo relation (Ziegler & Bender 1997), that cluster 
ellipticals form a homogeneous class of old objects. Pushing 
the observations to high redshift clusters allows to tighten 
the constraint, but only for big objects, not much smaller 
than 10^^ Mq. Besides, many ellipticals, presumably field 
objects, are predicted to form at lower redshift, z < 2. This 
is consistent with observations, which suggest that field el- 
lipticals are not a separate class of objects, but are younger 
on average (Bernard! et al. 1998; Franceschini et al. 1998; 
Abraham et al. 1998; Manenteau et al. 1998). 

The inversion of hierarchical order is in line with the 
general trend of galaxy formation: the star-formation his- 
tory at z < 1 is dominated by dwarf galaxies, while L* 
galaxies, both elliptical and spiral, seem already in place at 
z = 1 (see, e.g., Ellis 1998). However, a direct determina- 
tion of a 2/ — Mh relation from the age of stars in ellipticals 
is hampered by the well-known age-metallicity degeneracy. 
When trying to break this degeneracy, a possible age depen- 
dence, consistent with the one proposed here is claimed by 
many authors (Matteucci 1994; Bressan, Chiosi & Tantalo 
1996; Franceschini et al. 1998; Caldwell & Rose 1998; Fer- 
reras. Chariot & Silk 1998; Pahre, Djorgovski & de Carvalho 
1998). 



2.5 Predictions on galaxies 

Given the number of galactic halos nniMH ,p, Zsh) (Eq^ ^ 
HHi P and |lo|) , calculated at the shining time Zsh, the total 
number of galaxies at redshift z — with morphological type 
defined by pi < p < P2 (the limits can depend on Mh) is 
readily calculated by integrating the number density of halos 
in p and Zsh- 



Igal(MH) dMH = 



(12) 



dzsh / dp nH{MH,p, Zsh) ] dMH- 

Jpi 



It is supposed that the halo spin or the merging fraction 
are responsible for galaxy morphology. Both mechanisms are 
likely to have a role in determining whether a galaxy is going 
to be bulge-dominated: galaxy mergers give rise to "hot" 
galaxies, and the profiles of big ellipticals are consistent with 
the merging origin (see, e.g., Faber et al. 1997). On the other 
hand, low-spin systems can lead to large bulge-disc ratios, 
as sugg ested, e.g., by Mo, Mao & White (1997)[} 

Elliptical galaxies are assumed to be hosted either in 
low-spin (A < Xe) halos or in halos which have suffered a 
major merger (1/2 < / < /b). In both cases, the probability 
for a halo to host an elliptical increases with Mh/M^. As a 
consequence of the inversion of hierarchical order for galaxy 
formation, described in the previous subsection, small galac- 
tic halos shine later, when they are small with respect to M* 



t Galaxy morphology has probably a more complex origin: for 
instance, many lonticulars could come from spirals infalling into 
cluster; see, e.g., Ellis 1998 and references therein. 



at their shining time. Then, the fraction of elliptical galaxies 
is not fixed but increases with mass; this leads to a flatten- 
ing of the elliptical mass function with respect to the mass 
function of all the other halos. 

With a fixed p-threshold and a suitable tuning of the 
free parameters involved (see Section 2.7 for full details), 
it is possible to obtain a satisfactory prediction for the halo 
mass function of ellipticals, which is the main concern of the 
present paper. However, in Section 3.1 we will also test the 
predictions for the mass function of spiral halos, so as to give 
further support to the galaxy formation model presented 
here. 

The mass function of non-elliptical halos is steeper than 
that of spirals. This is in line with the general behaviour of 
galaxy formation models, which tend to predict a steep lumi- 
nosity function (see, e.g.. White 1993; Cole et al. 1994). This 
problem is usually solved by assuming a low efficiency of star 
formation for small halos. If spin is supposed to determine 
the galaxy type, it is possible to subtract small-mass, high- 
spin halos by assuming that they are not going to host bright 
galaxies but large low surface brightness (LSB) discs. Fol- 
lowing Dalcanton, Spergels & Summers (1997) and Jimenez 

et al. (1997)1^, the spin threshold As for the formation of 
non-bright galaxies is assumed to depend on mass: small- 
mass halos are more likely to host a non-bright galaxy. The 
spin threshold is parameterized as follows: 



log \s(Mh) = as log(Mff /lO^" M©) + log A^. 



(13) 



Then, a spiral halo is selected if Xe < A < Xs{Mh)- Con- 
sistency with observations is obtained for as ~ 0.4 (see Ta- 
ble 2), in rough agreement with Jimenez et al. (1997) who 
find As oc M^*, and Dalcanton et al. (1997) who report 
As oc Mjj^\ 

The predicted redshift evolution of the formation of 
bulge-dominated galaxies can be tested through the star for- 
mation history of elliptical galaxies given by Franceschini et 
al. (1998). It is possible to get a rough prediction of the 
contribution of ellipticals to the star formation rate of the 
universe by assuming that they form all their stars at the 
shining time of the halo, as defined in Section 2.4. Assuming 
that the visible mass of an elliptical, Mbui, is connected to its 
light through the relation Mbui/L = (M/L)o(I//L*e)'''"''"^ 
(where L*_b is the Schechter parameter for the E luminosity 
function), the bulge mass is related to the halo mass as: 



Mh 
ML 



/3bul//3B 



(14) 



The parameter /3b is the exponent of the Mh — L relation 
for ellipticals, and is defined in Section 3.1, Eq. ^; its as- 
sumed value is 0.75 (see also Table 1). Following paper I, the 
parameter /3bui is set to 1.25, while {M/L)q is set to 6.9/i 
(luminosities are in the B band). The star-formation rate is 
then calculated as: 



SFR(2) dt = 



dp dMnM^uiiMn) nHiMH,p,z) 

p^E Jo / 





dz 


) 


It 



(15) 



dt. 



■t According to Jimenez et al. (1997), high-spin halos do not even 
host LSB's, but remain dark. 
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Here the integral in p is performed over the interval which 
defines the elliptical morphology. 

Eq. ^ is only valid under the rather artificial hypothesis 
that all stars in an elliptical form in a very short time. This 
is a good assumption only at low redshift, when the age 
of the Universe is much larger than the typical duration of 
the starburst. A better estimate of the star-formation rate 
of ellipticals can be obtained by convolving Eq. |l^ with a 
curve describing an average star-formation history, given for 
instance by a truncated exponential with timescale r: 

SFR' {z)dt = J^ SFR(t,h) ^exp(^~ ^~^'^^^ ^ dteh (16) 

In this case, the shining time is identified with the first e- 
fold time, as the stabilization of the halo is supposed to 
be a result of the massive star formation, and then cannot 
precede it. In the following we will use for r a value of 1 Gyr, 
which is of order of the delay time of a 10^^ Mq halo. Note 
that this is probably larger than the star-formation timescale 
for giant ellipticals, which is likely to be smaller than 0.3 Gyr 
(Matteucci 1994). Of course, the star- formation timescale is 
physically related to the delay time tdoiay, and is likely to 
depend on halo mass. However, an accurate modeling of such 
a timescale is beyond the scope of the present paper. Then, 
the convolved star-formation history should be considered 
just as an indication of what can happen when relaxing the 
hypothesis of very fast burst. It is also noteworthy that, as 
expected, the convolution does not change appreciably the 
low-redshift cutoff of the star-formation history. 

2.6 Predictions on QSOs 

As mentioned in the Introduction, BHs strongly prefer ellip- 
tical morphologies, as their masses correlate with the mass of 
the host bulge. This correlation presents a sigiuficant scatter 
of 0.3-0.5 in decimal logarithm (see paper I). This scatter 
reveals the need of a "hidden variable" , able to modulate the 
efficiency of BH formation in halos, thus giving the wanted 
broadness in the BH-bulge relation. On the other hand, the 
mechanism responsible for galaxy morphology has an indi- 
rect infiuence on BH formation. As a reasonable working 
hypothesis, we assume that both efficiency of BH formation 
and galactic morphology are infiuenced by the same physical 
variable, spin or merging. 

The PDF of the morphological parameter must be such 
to reproduce, under a reasonable transformation, the high 
luminosity tail of the QSO luminosity function. Although 
merging is a physically motivated cause for stimulating the 
accretion onto a BH (it creates non-axisymmetric pertur- 
bations which help the gas to fall toward the center of the 
potential well), from the statistical point of view the merg- 
ing fraction / results unsuitable to this purpose: as shown 
in Fig. 1, the /-distribution is very fiat around the value 0.5, 
which corresponds to the merging of clumps of nearly equal 
mass, an event which is not asymptotically rare. Then, to 
shape this function into a steep power-law, so as to fit the 
QSO luminosity function, it is necessary to assume an ex- 
tremely steep dependence of the efficiency of BH formation 
on /. This is unphysical, as the fate of gas is not expected 
to be sensitive to infinitesimal variation of the ratio of merg- 
ing clumps. However, the merging fraction is not the only 
important physical quantity involved in a merging process. 



Other quantities, like the relative orientation of the spins of 
the merging clumps, or their impact parameter, are likely to 
be important in determining the final galactic morphology. 
Then, a more detailed modeling of the merging process is 
required to make it suitable for modulating BH formation, 
but this is beyond the scope of the present paper. 

The spin of the halo does not suffer from the problem 
discussed above, as low-spin halos are asymptotically rare. 
Halo spin has already been proposed as an important vari- 
able by Eisenstein & Loeb (1995b). On the other hand, what 
is physically relevant for the formation of a BH is the quan- 
tity of angular momentum that the gas is able to lose (see 
also Cavaliere & Vittorini 1998), and the internal distribu- 
tion of this angular momentum in the very center of the 
halo. Moreover, with a spin profile compatible with actual 
elliptical galaxies, it is very difficult to have the formation of 
a huge BH (De Felice, Yu & Zhou 1992). Then, the relevance 
of the global spin of the DM halo for BH formation is not ob- 
vious. But the BH-bulge relation shows that the final mass 
of the BH is related to some global property of the halo; 
then, it is not unreasonable to assume it to be related to an- 
other global property, as the total spin. The physical cause 
of this relation could be a dynamical feedback of the BH to 
the halo: a strong central mass infiuences through chaotic 
mixing the orbits which pass near the center, making them 
more axisymmetric (Merritt 1998). This could limit the mass 
of the BH, making it to depend on the total spin, which is 
relevant for axisymmetric systems. Alternatively, the BH- 
bulge relation could be due to the competition with star 
formation (Wang & Biermann 1998) or to the mechanical 
feedback of the BH on the protogalaxy (Silk & Rees 1998); 
the latter mechanism can give a steep dependence of the BH 
mass both on the bulge mass (M. oc M^^f) and on the halo 
spin (M. oc A~^) (Haehnelt, Natarajan & Rees 1998). In any 
case, a dependence of the BH mass on the total spin is ex- 
pected if this quantity infiuences the profile of the halo, and 
then its central density. Finally, as the halo spin is acquired 
from tidal torques given by the large-scale structure, a direct 
infiuence of it on the BH formation would explain the ob- 
servational evidence of alignment between radio loud AGNs, 
host elliptical galaxies and large-scale structure (West 1994). 

The mass of the BH changes rapidly in time during the 
bright QSO phase. However, as shown in paper I, the most 
important period, the longest and the brightest, is that in 
which the BH has acquired most of its final mass. It is then 
reasonable to consider only the final mass of the BH, M, , 
which is the result of the accretion onto the seed BH of 
the matter available in the reservoir. In the following, the 
instant at which the BH acquires its final mass will be called 
BH formation. Further accretion of mass in a non-bright or 
reactivated phase is assumed negligible. 

To construct a prediction for the QSO luminosity func- 
tion it is necessary to specify the efficiency of BH formation 
in DM halos: 

EH = M./Mh, (17) 

i.e. the amount of matter which ends up into the BH in units 
of the halo mass. This is assumed to depend on spin in the 
following way: 

EH = EHO i^J^J ■ (18) 
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The reference value Aq is set so as to be a 2a event when 
Mh = Aft, while eho and Oq are left as free parameters. 
Eq. ^ implies a linear scaling of M. with the halo mass 
Mh (neglecting the weak mass dependence induced by the 
spin). However, BH masses appear to scale with the bulge 
mass; this issue will be addressed in Section 5.3. 
The number of BHs formed at redshift z is: 



n,{M,,z)dM, dz = 

dM„ nHiMH,X{M.,M„),z] 



(19) 



dM, dz 

aq In 10 ■ 



To obtain a QSO luminosity from a BH mass, it is 
assumed that the BHs accrete mass at a fixed ratio of 
the Eddington rate, /ed = L/Led, where the Edding- 
ton luminosity is defined as Led = 4:TTGmpcM,/at ~ 
3.4 10* (M./Mq) Lq {mp is the proton mass, ar the Thomp- 
son cross section). Then, the BH mass grows to its final value 
during a timescale, or duty cycle time: 

iduty = eteD//ED, (20) 

where the Eddington time is defined a,s ts = M,c? / Ly,t> 
~ 4 10* yr and e is the efficiency of radiation of the QSO in 
units of Mc? , where M is the accreted mass. The luminosity 
of the QSO in a given e.m. band is related to the BH mass 
through: 



Lqso{M.) = /ed Led (M.)/Cs. 



(21) 



Here Cb is the bolometric correction appropriate for the 
e.m. band used. Observational evidence (see, e.g., Padovani 
1989; Wandel 1998) suggests that the efficiency of accretion 
is a function of the QSO luminosity, going from ~ 0.05 — 0.1 
for small AGNs to ~ 1 for bright QSOs. Following paper I, 
it is assumed that: 

/ r \ "ED 



/ed = 



\^1049erg/s 



(22) 



with the exponent qed set to 0.2. 

The luminosity function of QSOs is then: 



(23) 



M, 



dL 



QSO- 



nQSo(LQSo; z)dL = 

W.(A/.(Lqso), ^) (1 — QED)tduty 

Lqso 

It is noteworty that with the inversion of the hierarchical 
order for galaxy formation, the "inverted" evolution of QSOs 
is naturally obtained. 

The predicted mass function of dormant BHs at z = 
is obtained simply by integrating in redshift the number 
density of BHs given in Eq. M. 

The bivariate number density of BHs hosted in bulges 
of a given mass is obtained from Eq. ^ by transforming the 
halo mass and the spin into bulge (Eq. |lj) and BH (Eq. [l8| ) 
masses, and by integrating the resulting distribution in red- 
shift: 



n._bui(M.,Mbui)dM. dMbui = 

dz nH{M„(Mt,ui),HM.,M„),z] 

/3bu.a,lnl0 Mb.,M. 
© 1999 RAS, MNRAS 000, 
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This equation is valid only for halos which host elliptical 
galaxies. 

Finally, while it is reasonable to assume that the global 
properties of the galactic DM halo determine the amount of 
mass available to the BH for accretion, the actual amount of 
mass accreted could depend on some unpredictable details of 
the accretion process. This could infiuence the value of /ed, 
which would then become a random variable (Siemiginowska 
and Elvis 1997; however, this would infiuence the conclusions 
of paper I), or it could only infiuence the value of eho, which 
would then be modulated by a completely random process. 
These possibilities will be addressed elsewhere. 

2.7 Assumptions and parameters 

It is useful at this stage to list all the assumptions made in 
this Section: 

• The number of dark-matter halos dynamically formed 
at a given redshift, with a given mass and a given spin or 
merging fraction, is given by a set of "numerical recipes", 
mainly based on the extended PS formalism, which are 
known to fit N-body simulations. 

• Galactic halos are distinguished from those correspond- 
ing to galaxy groups and clusters according to a cooling cri- 
terion (which suppresses galaxy formation in large halos) 
and a dissipationless merging criterion (which suppresses 
galaxy formation at lower redshift). 

• Because of feedback, galaxy formation (or, in other 
words, the shining of the halo) is delayed with respect to 
the dynamical formation of the halo. The delay is larger for 
smaller halos. 

• The merging of galactic halos which have already shone 
is neglected. 

• Elliptical galaxies are hosted either in low-spin halos or 
in halos which have experienced a major merger. 

• High-spin halos do not host a bright galaxy. 

• The efficiency of BH formation in galactic DM halos is 
modulated by the same physical variable which is responsi- 
ble for the galactic morphology. 

• The QSO activity is close in time to the main burst of 
star formation for bulge stars. 

• QSOs shine only once, and acquire most of their mass 
in their bright phase. 

• The efficiency of accretion of the BHs, expressed in Ed- 
dington units, is a function of the BH mass; this is given in 
paper I. 

The free parameters of the model are listed in Table 2, 
together with the cosmological parameters which define the 
background cosmologies used. In the following Sections the 
parameters will be constrained by comparing the predictions 
of the model to many distributions inferred from observa- 
tions: the mass function of galactic halos at low redshift, 
divided in broad morphological classes, the star-formation 
history of elliptical galaxies, the luminosity function of (op- 
tical and obscured) QSOs and its evolution in redshift, the 
mass function of dormant BHs in nearby galaxies, and the 
scatter in the BH-bulge relation. For each cosmology, it is 
possible to find a set of acceptable parameters. 

The best fit parameters are obtained by means of a qual- 
itative comparison of the predicted and observed quantities. 
This is possible because each parameter infiuences mainly 
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some particular prediction, so that they can be fixed sepa- 
rately. In particular: 

• the a/ parameter (Eq. which determines the inver- 
sion of the hierarchical order for galaxy formation, is fixed by 
reproducing the correct slopes of the galaxy mass functions; 

• the normalization time tf (Eq. hll) is fixed by repro- 
ducing the high-redshift evolution of the QSO luminosity 
function; 

• the cutoff redshift zq (for Eq. ^|) is fixed by reproducing 
the cutoff of the star-formation history of elliptical galaxies; 

• A/cool (Eq. ^ is fixed by fitting the cutoffs of the galaxy 
mass functions; 

• \e and fs (limits of the p- integral in Eq. are fixed 
by reproducing the normalization for the mass function of 
ellipticals; 

• Ag (Eq. is fixed by reproducing the normalization 
for the mass function of spirals; 

• as (Eq. ^3|) is fixed by reproducing the slope of the 
spiral mass function; 

• the mass function of ellipticals in the merging case is 
sensitive to the value of the timescale A log b (Eq. |7|), which 
is set to 0.1, different values lead to unsatisfactory mass 
functions; 

• Oq (Eq. |l^) is tuned to obtain a good shape for the 
QSO luminosity function; 

• £ho (Eq. |l8| ) is obtained by fitting the mass function of 
dormant BHs; 

• e (Eq. is obtained by reproducing the normalization 
of the QSO luminosity function once eho is fixed. 



3 GALAXIES 

3.1 The halo mass function of galaxies 

Galaxies are embedded in DM halos, the mass of which 
is not directly observable; the statistical quantity which 
can be observed is the luminosity function. However, at 
fixed morphological type galactic halos are known to fol- 
low some kind of universal profile (Salucci & Persic 1998a). 
This is theoretically confirmed by the existence of a uni- 
versal profile for DM halos (Navarro, Frenk & White 1995; 
Moore et al. 1998), even though observational and theoret- 
ical halos differ in some details. This fact implies that halo 
masses and galaxy luminosities are strictly related through 
a morphology-dependent Mh/L relation, which is in general 
a function of luminosity. Once this relation is known, it is 
possible to determine the mass function of objects from their 
luminosity function, and to compare them to the predictions 
of the galaxy formation model of Section 2. 

We divide the morphological types for bright galaxies 
into two broad categories, one of early types (E and SO, 
briefly E), and one of late types (Sa to I, briefly S). The 
luminosity functions are parametrized by means of the usual 
Schechter formula: 

ML)dL = <j},^(L/L,,y°'^ exp{-L/L,^)dL/L,„ (25) 

where the index i is E or S. Type-dependent luminosity 
functions for E and S galaxies are given for instance by 
Efstathiou, Ellis & Peterson (1988), Loveday et al. (1992; 
Stromlo/APM survey), Marzke et al. (1994; CfAH~2 sur- 
vey), Heyl et al. (1997; the AUTOFIB survey), Marzke et 





E 


S 


</>.(Mpc-3h3) 


3.8 10-3 


1.1 10-2 


a 


1.0 


1.2 


Mg{mag — 5 log h) 


-19.8 


-19.8 


Mh*(Mq), EdS, h = 0.5 


1.45 10^2 


1.45 10^2 


Mh*{Mq), Lambda, h = 0.7 


0.67 10i2 


0.67 10^2 


Mh*(Mq), open, h = 0.7 


0.67 1012 


0.67 1012 


P 


0.75 


0.56 



Table 1. Parameters for the galaxy luminosity functions and 
Mff/L relations {h = 1 unless otherwise stated). 

al. (1998; the SSRS2 survey), Marinoni et al. (1998, NOG 
sample). There is a broad agreement on the values of the 
various parameters, but different authors disagree in some 
important details. In particular, the E luminosity function 
may be flatter than the S one, as suggested for instance by 
the early work of Efstathiou, Ellis & Peterson (1988), or by 
the AUTOFIB survey, but not confirmed by the CfAl-|-2 or 
SSRS2 surveys. A recent reinvestigation based on a large lo- 
cal sample of galaxies (Marinoni et al. 1998) has confirmed 
that the slope of the luminosity function steepens gradually 
from E to Sm/I galaxies, with the exception of SO galax- 
ies, whose luminosity function is as steep as that of Sc-Sd 
galaxies. It is then apparent that the true slope depends 
sensitively on the definition of the sample used. 

In this work, values of 1 and 1.2 will be used for the 
slopes of the E and S luminosity functions; these values are 
consistent with most determinations. Table 1 shows the val- 
ues adopted for the parameters of the E and S luminosity 
functions, which are roughly consistent with all the lumi- 
nosity functions listed above. 

Following Salucci & Persic (1998a), the mass- luminosity 
relations are assumed to be of the kind: 

Mh ^ Mh,{L/L„f', (26) 

where again i is E or S. Spiral galaxies have a /3s parameter 
of 0.56 and an = 2.4 10^^ M0/i-i(n(zsh)/no)~^/^(l + 
ZBh)~^ (Persic, Salucci & Stel, 1996), where z^h is the average 
shining redshif of the halo. We take values 2.3, 2.5 and 3.0 
for the EdS, Lambda and open models; the results, reported 
in table 1, do not depend much on the exact value of this 
parameter. The Pe and Afjf^, parameters, relative to early- 
type galaxies, are much harder to obtain, as the evidence of 
DM is detected only in the outer regions (see, e.g., Danziger 
1998). The limited evidence available suggests a /3b param- 
eter not so different from /3s, and most likely smaller than 
one (Salucci & Persic 1998a). The value 0.75 will then be 
used in the following. The Afjf, parameter is set equal to 
ML- 

Fig. 3 shows the comparison of the predicted to the 
"observed" mass functions in the cases of EdS and Lambda 
Universes (the open case is very similar to the Lambda one, 
and is not shown). Table 2 gives the best-fit values of the 
parameters used, for the different cosmological models. Al- 
though ellipticals are the main concern of the present paper, 
we show also the results for spirals in order to give further 
support to our model and to our criterion for identifying the 
morphological type. 

In the EdS case, the mass functions are reproduced with 
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Figure 3. Mass functions of galactic halos at z = 0. 



roughly the correct slope and normalization. The elliptical 
mass function obtained with the merging fraction is steeper 
than the spin one. The exponential cutoffs are not perfectly 
reproduced, and in the spin case ellipticals do not manage 
to outnumber spirals at large mass; this is due to the fact 
that no inversion of hierarchical order is present at large 
masses (Fig. 2). We have chosen to tune Mcooi, so as to best 
reproduce the mass function of ellipticals. 

In the Lambda case, the overall normalization is slightly 
underestimated. However, the discrepancy is within the er- 
ror in the normalization of the observational mass func- 
tions, which is influenced by the uncertainty in the model- 
dependent quantity M§,, and by the Hubble constant. The 
predicted slopes for the mass functions tend to be slightly 



steeper than in the EdS the fraction Mh/M,{z) 

varies more slowly in a non-critical Universe; this is cor- 
rected by increasing the parameter a/ to 2.6. Again, as el- 
lipticals are more relevant in this context, we have chosen to 
tune the Xe and fc parameters so as to reproduce at best 
the correct number of ellipticals, at the expenses of spirals. 

The mass function of spiral halos (in the spin case) de- 
pends on the subtraction of LSB halos, as explained in Sec- 
tion 2.5. We have verified that the predicted mass function 
for LSB is consistent with the luminosity function given by 
Sprayberry et al. (1997) and the M/L ratio suggested by 
Salucci & Persic (1998b). These aspects of galaxy formation 
will be discussed elsewhere. 

It appears that both spin and merging give mass func- 
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Figure 4. Star formation history of elliptical galaxies, with an 
assumed timescale r of 1 Gyr (Eq. [l6[ ) , compared to the result of 
Franceschini et al. (1998) (shaded region), and, for reference, to 
the data of Madau (1997) (squares; corrected for dust extinction 
as described in the reference) and Steidel et al. (1998) (stars). 



Table 2. Values of the free parameters for the cosmological mod- 
els considered. 

tions which are roughly consistent with the ones inferred 
from observations. This implies that at this stage we do not 
need to choose between the two mechanisms. This conclu- 
sion will change in Section 4. 

3.2 The star-formation history of elliptical 
galaxies 

Fig. 4 shows for the EdS case the cosmological star forma- 
tion rate of ellipticals as inferred by the analytical model, 
assuming that bulge stars form with a truncated exponential 
history with timescale 1 Gyr (Eq. ^ see Section 2.5). It is 
worth recalling that this curve is considered as just indica- 
tive of what happens when the hypothesis of very fast burst 
is relaxed. The prediction is compared to the estimate of 
Franceschini et al. (1998), relative to ellipticals in the Hub- 
ble deep field, and for reference to the estimates of the global 
star-formation history by Madau (1997) (corrected for dust 
extinction as discussed in the reference) and Steidel et al. 
(1998). 

There is an inconsistency between the Franceschini et 
al. (1998) curve and the bulge mass function given in paper 
I: the latter gives a mass density of 6.3 10^ Mq/Mpc~^ in 
bulge stars, while the former gives from 1.2 10* to 2.5 10* in 
the same units. The model reproduces a halo mass function 
which is consistent with the bulge mass function given in pa- 
per I, and then underestimates the Franceschini et al. curve 
by a factor from 2 to 4. Taking this inconsistency into ac- 
count, the predicted star-formation history is in agreement 
with that of Franceschini et al. (1998): it shows a decline at 
z < 2 and a very broad peak at z; ~ 3, in rough agreement 



with the flatness of the star-formation history suggested by 
Franceschini et al. (1998) and Steidel et al. (1998) (see also 
Pascarelle, Lanzetta & Fernandez-Soto 1998), even though 
the normalization is not recovered. The star formation his- 
tory of ellipticals will be addressed in more detail in a forth- 
coming paper (Silva et al. 1999). 

On the other hand, the star-formation history of all 
stars at < z < 1 is not reproduced; star formation is dom- 
inated at low redshift by stars in spiral discs and in dwarf 
galaxies, i.e. objects which are not related to the QSO phe- 
nomenon. Note that this conclusion is at variance with Boyle 
& Terlevich (1998). 

3.3 Relationship with z ~ 3 galaxies 

In the last years, an increasingly large sample of galaxies at 
2 ~ 3 has been observed; such galaxies are actively star- 
forming galaxies, found as "UV dropouts" in deep flelds 
(see, e.g., Steidel et al. 1998). These are often referred to 
as Lyman break galaxies (hereafter LBG). The clustering 
properties of such galaxies are consistent with a scenario 
in which each galaxy is associated with a single halo of 
~ 10^^ A/0 (Adelberger et al. 1998); direct dynamical es- 
timates of halo masses seem to indicate smaller values, but 
the results are still too uncertain to draw firm conclusions 
(Pettini et al. 1998). It is interesting to note that the clus- 
tering of QSOs, which evolves slowly with redshift, is more 
similar to that of LBGs than to that of general galaxies 
(La Franca, Andreani & Cristiani 1998; Magliocchetti et al. 
1998). The abundance of LBGs is 6.4 10"^/i^Mpc"^ for the 
EdS case, and 1.7 10"^/i^Afpc"^ for an open model (Stei- 
del et al. 1998). Such galaxies have been interpreted as the 
ancestors of big ellipticals (Governato et al. 1998). 

According to Fig. 2, LBGs would correspond to galaxies 
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forming in halos of ^ 10^^ Mq, in broad agreement with the 
estimate given above. As a substantial fraction of L* galaxies 
are forming at that redshift, the predicted abundance of such 
halos is of the order of $* for ellipticals; if it is assumed that 
only bulge-dominated galaxies are actively forming stars at 
those redshifts, then the abundance of halos able to host 
LBGs is broadly compatible with the values given above (see 
table 1). This implies that LBGs do not sample the whole 
population of DM halos, but only those whose characteristics 
(spin or merging in this context) are such to cause strong 
star formation, while a majority of halos will host proto- 
spirals which are not visible at that redshift. 

In the present context, as long as each LEG is associ- 
ated to a single halo, these objects are going to survive to the 
present epoch as big ellipticals; further mergings should be 
negligible in terms of mass and star formation. Then, the sce- 
nario presented here is not in complete agreement with Gov- 
ernato et al. (1998), KaufFmann, Nusser & Steinmetz (1997) 
and Baugh et al. (1998), who predict that Lyman break 
galaxies are just pieces of big ellipticals, which must merge 
subsequently. The disagreement is weakened if the one-to- 
one correspondence of halos and galaxies is relaxed; in this 
case the merging of LBGs contained in a single halo does not 
imply any halo merging. An interesting example is given by 
the observation of the distant radio galaxy 1138-262 (Pen- 
tericci et al. 1997, 1998): while optical observations reveal a 
large amount of substructure, corresponding to many (about 
ten) Lyman-break objects, observations in the (rest frame) 
NIR and Lyman-a emission show a much more coherent 
structure, with a velocity dispersion of ~ 300 km/s. This 
is consistent with a young giant elliptical, with some knots 
of star formation visible in the optical as separated enti- 
ties. Multi-band observations of many elliptical galaxies are 
needed to assess the dynamical state of high-redshift galax- 
ies. 



4 QSOS 

4.1 Observational properties 

The optical luminosity function of QSOs is well known (see, 
e.g., Boyle, Shanks & Peterson 1988). The typical luminos- 
ity of QSOs evolves rapidly with redshift: at smaller redshift, 
Si 3, the evolution is similar to a pure luminosity evolu- 
tion with Lt{z) oc (1 + 2;)^'^ for z < 2, and Lt{z) ~ const at 
2 < 2 < 3, while at z > 3 QSOs start to decrease in number 
(see, e.g., Osmer 1998, Shaver et al. 1998). At smaller red- 
shift, z < 1, the evolution is not purely in luminosity, as the 
luminosity function tends to flatten (La Franca & Cristiani 
1997). 

To describe the shape and evolution of the QSO lu- 
minosity function, the parameterization proposed by Pei 
(1995) has been used (for a spectral index a = —0.5). Lumi- 
nosities are given in the B band; following paper I and Elvis 
et al. (1994), the bolometric correction has been set to 13. 

We assume the existence of a whole population of ob- 
scured AGNs, which contribute to the cosmological back- 
ground in the hard X-rays (Setti & Woltjer 1989; Celotti et 
al. 1995; Comastri et al. 1995). Following Comastri et al. 
(1995), the abundance of such objects is estimated starting 
from the soft-X-ray luminosity functions of QSOs (Boyle 



et al. 1993). It is assumed that for each observed soft-X- 
ray QSOs there are 5.4 ones which are heavily obscured. 
The applied bolometric correction is 25 (see again paper I 
and Elvis et al. 1994). It results that obscured QSOs dom- 
inate the mass function of dormant BHs at small masses, 
M, < 10* Mq, but do not contribute at larger masses. 

Figs. 5 and 6 show (for the EdS and Lambda models) 
the optical luminosity function, the estimated contribution 
of obscured objects and their sum; both are given in terms 
of the bolometric luminosity. Although the calculations have 
been performed using the analytical parametrizations, the 
data points from Pei (1995) and Boyle et al. (1993) are 
shown in the flgure. At z — 4.4, the newer data points 
of Kennefick, Djorgovski and Meylan (1996) are reported 
in place of the Pei data ai z — 4. Obscured objects con- 
tribute significantly only in the low-luminosity end, while 
bright QSOs are almost unobscured. The contribution of 
obscured QSOs is not considered at z > 3, as the X-ray 
luminosity function is not measured there. It is interesting 
to see that the composite luminosity function appears al- 
most featureless and quite steep at all luminosities. In this 
case, the knee of the optical luminosity function would be 
interpreted as an effect of the onset of obscuration. It must 
be noted that both luminosity functions are still uncertain, 
which makes the complete luminosity function still specula- 
tive. Nonetheless, the qualitative trends of the dominance of 
obscured objects only at small luminosities and the steep- 
ening of the complete luminosity function should be robust. 

The evolution of the luminosity function is quantified, 
following Pei (1995), through the index (I/boi) 8- function 
of redshift. The exponent 2 is such to give more weight to 
the objects near the knee of the luminosity function, which 
is most robust; in this way the evolution index is not sensi- 
tive to errors at both ends of the luminosity function. This 
index is not sensitive to the contribution of obscured ob- 
jects, which is then neglected. The evolution index (iboi) is 
shown in Fig. 7. 

4.2 Comparison of predictions with observations 

Figs. 5, 6 and 7 show the comparison of the predicted and ob- 
served QSO luminosity functions and evolution index (Lboi) > 
for the EdS and Lambda models. The free parameters in- 
volved are e, a, and eho (Section 2.6), together with tf and 
a/ (Section 2.4); the best-fit values are again given in Table 
2. The parameters have been chosen so as to fit the com- 
plete luminosity function; as the predicted curve is almost 
featureless (it is approximately a rescaling of the low-spin 
tail of the lognormal spin PDF), it is much easier to fit sat- 
isfactorily the complete luminosity function rather than the 
optical one only; anyway, a moderately good fit of the opti- 
cal curve can be obtained by increasing the value of the a, 
parameter. 

The agreement between model and data is overall very 
good. At smaller redshift, z < 1, the number of QSOs 
is slightly underestimated, especially at small luminosities. 
This modest disagreement may imply that recurrency of 
AGN activity is present at small redshift. At high redshift 
the model predicts a luminosity function which is fiatter 
than the one extrapolated from lower redshift, while lower- 
luminosity activity is severely suppressed. This is consistent 
with the currently available upper limits of Kennefick et al. 
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Figure 5. QSO luminosity functions at different redshift. EdS case. The optical luminosity function is taken from Pei (1995), the 
contribution of obscured objects is based on the X-ray luminosity function of Boyle et al. (1993). The curve named "ALL" gives the 
total contribution of optical and obscured objects, and is based on the parametrizations of the optical and X-ray luminosity functions 
given by the already mentioned authors. At 2; = 4.4 the newer data points of Kennefick et al. (1996) (denoted by stars) are reported. 



(1996), although the EdS prediction tends to overpredict the 
high- luminosity tail. The open model, not shown for brevity, 
tends to predict a larger number of high redshift QSOs, un- 
less the delay time tf is more than 0.6 Gyr. 

The parameter a, is set to 1.8, 1.3 and 1.7 in the EdS, 
Lambda and open cases, confirming that spin should play a 
major role in BH formation, but it cannot be as large as 5, 
the value suggested by Haehnelt et al. (1997). The efficiency 
of radiation turns out to be equal to the canonical value 0.1; 
this is a good consistency test for the model. 



5 DORMANT BLACK HOLES 

5.1 The mass function of dormant BHs 

The determination of the mass function of dormant BHs in 
ordinary galaxies at z = heis been addressed in full detail 
in paper I, and briefly described in the Introduction. Here 
we present a comparison of the mass functions as obtained 
from the QSO luminosity function and the radio luminosity 
function of elliptical cores. The mass function of dormant 
BH masses is obtained from the QSO luminosity function (it 
is called AMF, as in paper I) under the same assumptions 
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Figure 6. The same as Fig. 5 for the Lambda case. 



described in Section 2.6, i.e. of accretion at /ed (Eq. |22[ ) 
times the Eddington limit for a time tduty (Eq. Eot): 



n,{M,)dM, = 



In 10 



(1 — Oed) tduty 



(27) 





dz 








IE 



nQso{LQso{M.)- z) dlogM. 



Obscured QSOs have been included by using the complete 
luminosity function for nQso{L; z). 

The radio-based mass function (called RMF) was ob- 
tained, as in paper I, by transforming the radio core lumi- 
nosity function of ellipticals (by Sadler et al. 1988, corrected 
to subtract non-core emission as in paper I) . The mass func- 



tion obtained in paper I (OMF in that paper) from the bulge 
mass function is not shown, as it is trivially satisfied when 
both the mass function of elliptical galaxies and the effi- 
ciency of BH formation (Section 5.2) are correctly repro- 
duced. 

Fig. 8 shows the comparison between the predicted BH 
mass function and the AMF and RMF. Again, the EdS and 
Lambda models are shown. The agreement is again very 
good. This result is trivial once the QSO luminosity function 
is correctly reproduced at any redshift, and once the consis- 
tency of the AMF and RMF has been assessed in paper I. 
However, fitting the AMF and RMF allows one to determine 
the parameter eho separately from e, as this quantity does 
not depend on the light actually emitted by the QSO. 
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Figure 7. Evolution index (L^ .) for tiio QSO population. Data arc taken from Pei (1995). Left: EdS case; right: Lambda case. 
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Figure 8. Mass functions of dormant BHs. The distribution n,{M,) is in dlogA/, 



Fig. 8 shows also the predicted contribution to the BH 
mass function of elhptical galaxies alone (defined through 
the spin threshold), which is the quantity to be compared 
to the RMF; the contributions with and without ellipticals 
differ where the RMF is not defined. As in paper I (see 
their Fig. 5), ellipticals give the dominant contribution to 
the mass function for M, ^ 10* Mq , while spirals dominate 
at smaller BH masses. This gives further support to the 
threshold criterion used to separate ellipticals from spirals. 



5.2 The M./Mbui ratio. 

A more interesting test would rely on the prediction of the 
joint number density of BH and bulge masses, given by 
Eq. This bivariate distribution gives full information on 
the correlation between BH and bulge masses, and could 
be compared to direct estimates of BH masses in nearby 
early- type galaxies. However, the available samples of galax- 
ies with known BH masses are just compilations of galax- 
ies for which suitable observations were obtainable, and the 
dynamical measures of BH masses are still affected by sys- 
tematic uncertainties connected to the kinematical models 
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used (Magorrian et al. 1998; van der Marel 1998; Ho 1998). 
As a consequence, it is not possible at present to obtain a 
reliable observational estimate of the bivariate M. — Afbui 
distribution. 

A more robust test relies on predicting the PDF of the 
quantity: 

7? = M./Mbui. (28) 

The ratio R is computed from the joint number density of 
BH and bulge masses (Eq. ^) by integrating it over one vari- 
able. The PDF of R has been estimated by Magorrian et al. 
(1998), on the basis of their sample of BH masses; they fit 
their PDF with many functions, among which a truncated 
decreasing power-law and a lognormal turn out to be ac- 
ceptable fits. Paper I gives a different estimate of this PDF: 
assuming a lognormal shape, the distribution must be such 
to obtain a good BH mass function from the mass function 
of galaxy bulges (the OMF), consistent with the AMF and 
RMF. The obtained values for the width and mean of the 
lognormal are 0.3 (in decimal logarithm) and 10~^'^; the 
width is slightly smaller than the one recovered by Magor- 
rian et al. (1998) (0.5 in decimal logarithm), implying that 
random errors in the estimates of BH masses are not a main 
source of the scatter in the M, — Mbui relation. As already 
mentioned in the Introduction, the average obtained in pa- 
per I is smaller than the Magorrian et al.'s value by a factor 
of ~ 2, in agreement with van der Marel (1998) and Ho 
(1998), and this is commented in full detail in paper I. 

The BH mass has been supposed to scale with the 
halo mass, which implies a BH-bulge relation with slope 
Ps/Pbui ~ 0.6. As a consequence, the R PDF depends much 
on the range of bulge masses over which it is averaged. The 
estimate of the R PDF given in paper I is sensitive only to 
the large-mass part of the mass function, as it is designed 
to spread the sharp cutoff of the bulge mass function into 
the milder cutoff of the BH mass function. Besides, the BH- 
bulge correlation is best tested for big ellipticals, while at 
lower bulge masses the measures are very uncertain and spi- 
ral bulges come into play. Then, the predicted R PDF is 
calculated only for ellipticals corresponding to luminosities 
larger than L,. Fig. 9 shows the comparison of the predicted 
R PDF and the one estimated in paper I. In both the EdS 
and the Lambda cases, the peak is reproduced at the cor- 
rect position, while the distribution is skewed toward large 
R values and truncated at low R values by the spin thresh- 
old. In case of poor and noisy data, this distribution would 
be easily fit by the lognormal given in paper I. Notably, at 
large R values the slope of the PDF is consistent with that 
of the truncated power-law of Magorrian et al. (1998), which 
is also shown in Fig. 9. It must be stressed that at this stage 
the R PDF is reproduced without tuning any parameter. 
The parameter which mostly influences the shape of this 
curve is a^, which is already constrained by the QSO lumi- 
nosity function. In other words, a connection is estabilished 
between the shape of the QSO luminosity function and the 
shape of the R PDF. 

As mentioned above, the mass of formed BHs has been 
supposed to scale linearly with the halo mass, and this im- 
plies a BH-bulge relation with slope ~ 0.6. On the other 
hand, the data of Magorrian et al. (1998) seem to suggest a 
linear or steeper relation between BH and bulge mass, while 
van der Marel (1998) suggests that the relation could be 
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Figure 9. PDFs of the R ratio defined in Eq. ^ The thick 
continuous lognormal curve is obtained from paper I, while the 
thin continuous line is the power-law fit of Magorrian et al. (1998). 
The EdS and Lambda curves are shown. 

shallower than linear, in line with our prediction. To force 
a linear BH-bulge relation, one can simply assume that the 
efficiency of BH formation scales with the bulge mass. This 
is implemented by multiplying the right-hand-side of Eq. |l^ 
by a term {Mh/{W^^ Mq))""'^, so that BH masses scale 
as M^", then setting an = Phui/PE- With this hypothesis, 
it is possible to reproduce satisfactorily the QSO luminos- 
ity function and dormant BH mass function; in the EdS 
model this is done by lowering the Uq parameter to 1 and 
the £ho parameter to 10~^'^. However, the a, parameter 
being lowered, the resulting R PDF is significantly narrower 
than those shown in Fig. 9, and then not compatible with 
the observational evidence. 

Then, in the present framework a linear scaling of M. 
with A/bui is not consistent with the data, unless a sig- 
nificant part of the scatter in the M. — A/bui correlation 
is due to observational errors, or unless a further mecha- 
nism is able to increase the scatter in the BH-bulge relation, 
without changing the QSO luminosity function; dissipation- 
less galaxy mergers which take place after the QSO phase 
could provide such a mechanism. An even steeper relation, 
M, oc Af^^j^, is suggested by Silk & Rees (1998); such a steep 
dependence would be hard to reconcile with the framework 
presented here. 



6 SUMMARY AND CONCLUSIONS 

This paper describes an analytical model which addresses for 
the first time the joint formation of galaxies and QSOs. The 
model is able to predict the halo mass function of galaxies 
of different broad morphological types, the star-formation 
history of elliptical galaxies, the QSO luminosity function 
and evolution, the mass function of dormant BHs in nearby 
galaxies, and the correlation of these with the mass of host 
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bulges. The model has been successfully compared to avail- 
able observational data, so as to constrain its free parameters 
(listed in Table 2). An acceptable fit is found for each of the 
cosmological models considered (EdS, Lambda, open). 
The following conclusions can be drawn: 

• Consistency with observation is obtained if the hierar- 
chical order is inverted when considering the shining epoch 
of galactic halos: large galactic halos experience massive 
star formation and QSO activity before smaller ones. This 
is achieved by delaying the shining of small objects. The 
inversion of hierarchical order is consistent both with the 
evolution of QSOs and with evidences based on stellar pop- 
ulations of elliptical galaxies. 

• There is a clear need of a "new variable" , which modu- 
lates the efficiency of BH formation for halos of a given mass. 
As QSOs prefer early-type morphologies, it is assumed that 
the "new variable" is the same as the one responsible for 
galaxy morphology. 

• The spin of DM halos is a good candidate as "new vari- 
able". The merging fraction appears a good candidate for 
determining galaxy morphology, but not for modulating the 
efficiency of BH formation, as the merging of halos of similar 
size is not asymptotically rare as bright QSOs are. A more 
detailed description of mergers could solve this problem. 

• The inversion of the hierarchical order of galaxy for- 
mation, together with a threshold criterion for predicting 
morphological types (based either on spin or merging) pro- 
duces mass functions for galactic halos which have the cor- 
rect morphology-dependent slopes. 

• A connection is estabilished between the slope of the 
QSO luminosity function and the distribution of the ratio 
M./A/bui. 



the process, and its agreement with many different pieces of 
observational evidence is encouraging. 

Our analysis shows that a delay of the "shining phase" 
of halos, when star formation and QSO activity occur, makes 
it possible to explain the main statistical features of the joint 
QSO-galaxy formation. This delay leads to an inversion of 
the hierarchical order for galaxy formation. This highlights 
the potential impact that QSOs can have in galaxy forma- 
tion. Their importance relies on the ever growing evidence 
that QSO activity is intimately related to galaxy forma- 
tion; the BH-bulge relation is a striking demonstration of 
such a relationship. The inversion of hierarchical order is 
not in contradiction with standard hierarchical CDM mod- 
els, as it is related not to the dynamical formation but to 
the shining phase of DM halos, and it is supposed to be 
caused by feedback mechanisms. The inversion of hierarchi- 
cal order leads to the prediction that smaller galaxies are 
made up of younger stars. This prediction can be tested by 
observations which are able to break the well-known age- 
metallicity degeneracy which affects old stellar populations; 
the evidence already available is consistent with it (Mat- 
teucci 1994; Bressan, Chiosi & Tantalo 1996; Franceschini 
et al. 1998; CaldweU & Rose 1998; Ferreras, Chariot & Silk 
1998; Pahre, Djorgovski & de Carvalho 1998). 
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A caveat on the role of spin is necessary: as only the sta- 
tistical properties of spin have been used in the model, the 
results presented here do not give direct evidence that spin 
is in play in determining the luminosity of QSOs, but reveal 
that any physical variable which has the same behaviour as 
the spin, i.e. has a PDF nearly lognormal in shape slowly 
changing with Mn/Mt,, is a viable variable for both deter- 
mining the morphological type and modulating the mass of 
the BHs. 

This paper gives further support to the favoured sce- 
nario of paper I, in which the most luminous QSOs (asso- 
ciated to BHs with masses larger than 10* Mq) are hosted 
in elliptical galaxies, shine only once for an Eddington time 
and at the Eddington limit, and are hardly obscured, while 
fainter AGNs are associated to spirals, shine at a fraction 
of the Eddington limit, may be significantly reactivated and 
are often heavily obscured. 

The analytical model presented here is based on a num- 
ber of reasonable but simplified assumptions. A more de- 
tailed description of DM halos would require the use of 
semi-analytic merging trees, or of large N-body simulations. 
Merging of galactic halos after the QSO phase is neglected 
by construction. Spin and merging are treated as alternative 
quantities, while they are likely to have both a role in shap- 
ing galaxies and triggering BH formation. The reactivation 
of existing QSOs, important to reproduce the low-level activ- 
ity at low redshift, is neglected. Nonetheless, the analytical 
model is supposed to catch the most important elements in 
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